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Abstract. A complete description of the nucleon structure in terms of generalized parton distributions 
(GPDs) at twist 2 level requires the measurement/computation of the eight functions H, E, H, E, Ht, 

Et, Ht and Et, all depending on the three variables x, £ and t. In this talk, we present and discuss our 
first steps in the framework of lattice QCD towards this enormous task. Dynamical lattice QCD results 
for the lowest three Mellin moments of the helicity dependent and independent GPDs are shown in terms 
of their corresponding generalized form factors. Implications for the transverse coordinate space structure 
of the nucleon as well as the orbital angular momentum (OAM) contribution of quarks to the nucleon spin 
are discussed in some detail. 

PACS. 12.38.Gc Lattice QCD calculations - 14.20.Dh Protons and neutrons 



1 Introduction to GPDs 

The generalized parton distributions (GPDs) of quarks in 
the nucleon we would like to compute are defined through 
the parameterization of off-forward nucleon matrix ele- 
ments 

(i y A ,| j g e ^ ( _A n)7 ^ ( A n)|P)A) 

= U(P', A') [j»H(x, t) + l -^±^E{x, f, t)J U(P, A),(l) 

for the helicity independent respectively vector case, and 

(P', A'| J g e ^(-\)7 57 ^(\)|P,A) 

= U(P', A') (j 5 ^H(x, f , t) + -^E(x, t t)j U{P, A),(2) 

for the helicity dependent /axial vector case Here 
we dropped for simplicity the dependence on the resolu- 
tion scale Q 2 as well as the gauge links U rendering the 
bilocal operators gauge invariant. The definition of mo- 
menta and helicity labels can be inferred from Fig.ljl"]). 
and in addition we note that A = P' — P, t = A 2 , while 
n is a light-cone vector with n± = 0, so that the longitu- 
dinal momentum transfer is given by £ = — n ■ A/2. The 
proper definition of the tensor/quark helicity flip GPDs 
can be found in 0]. Some of the excitement about GPDs 

* Based on talk given by Ph.H. 
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Fig. 1. Graphical representation of GPDs as part of a scatter- 
ing amplitude 



is due to their unifying nature, which can be read off right 
away form the definitions Eas. lll2fl . There is the forward 
limit, A — > 0, in which E and E do not contribute and 
the above definitions reduce exactly to the well known 
parton distributions, H(x,£ = 0, t = 0) = q(x) = fi(x), 
H{x,^ = 0, t = 0) = Aq(x) = gi(x). Secondly, by in- 
tegrating over x, we recover the defining equations for 
the the Dirac and Pauli form factors Fi (t) , F-i (t) as well 
as the axial and pseudo-scalar form factors gA(t), gp(t) 
(all to be taken with an appropriate flavor combination), 
$dxH{x,£,t) = Fi(t), J dxE(x,£,t) = F 2 {t) etc. Sim- 
ilar statements hold for the tensor GPDs as well. It is 
known for a long time that electromagnetic form factors 
encode information about the distribution of charge in the 
nucleon in terms of e.g. charge radii, upon Fourier trans- 
formation in the momentum transfer A to coordinate or 
impact parameter space, b. The obvious question is if and 
how this translates to GPDs. Due to the Heisenberg un- 
certainty relation, it seems at first sight that the GPDs 
as functions of the quark longitudinal momentum fraction 
x and A will not allow for a probability interpretation, 
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since longitudinal momentum and position in longitudi- 
nal direction cannot be determined simultaneously. As M. 
Burkardt has pinpointed [5], it turns out that just by in- 
tegrating over longitudinal position, that is setting £ = 0, 
ones gets functions of x and b± which have a perfect prob- 
ability density interpretation. To be precise, it is shown 
that 

q(x,b±) = i?(x,0, b±) 

d 2 A_ L e- ibxA - L H(x,£ = 0,t = -A 2 ± ) (3) 

is positive definite and interpreted as helicity independent 
quark probability density in b± (and x), even in a relativis- 
tic framework. Similar statements hold true for the impact 
parameter dependent helicity and tensor/transverse spin 
distributions, Aq(x,b±) and Sq(x,b±). It is quite illumi- 
nating to rewrite q(x, b±) using a fock state decomposition 
and nucleon wave functions W n c (overlap representation) (21 

EU, 



procedure, we get the following conservation laws (the first 
one also frequently called nucleon spin sum rule)|llj. 



9(^) = EE( 47r )" _1 / 



j J dxj d A 



x5 2 \b±-(l- x)r ±a + ^ Xjr±j 
x \& n>c (xi,...;r±i,...)\ 2 , 



(4) 



2=^ 



J 9 = - 2 (A 20 (0) 



B 20 (0)) 



g+g 



l = A q + 9 (0), = B q + 9 (0). 



(5) 
(6) 



where n is the number of partons in a given Fock state and 
the sums are carried out over all contributing n-parton 
configurations; a denotes the active parton. Eq.Q) shows 
clearly the anticipated positivity of the impact parame- 
ter distribution. Furthermore it becomes clear that de- 
tailed knowledge of the GPDs will certainly lead to an im- 
proved understanding of the all-important nucleon wave 
function. The first two delta functions in Eq.Q ensure 
conservation of longitudinal and transverse momentum re- 
spectively that the center of momentum (CM) is R± = 
y~]j—i Xj r ±j = 0. The last delta function shows on the one 
hand that b± is the distance of the active parton to the 
CM, b± = r± a — Y^j x j r ±j = r ±a- On the other hand, 
written as in Eq.Q, we see explicitly what happens if x 
is close to one (while Xj^ a — > 0). Taking into account that 
the distribution will finally vanish in the limit x — » 1, we 
observe that the normalized distribution becomes strongly 
peaked at the origin, lim x ^>iq(x,b±)/q(x) oc S 2 (b±). In 
momentum space, we get the remarkable prediction that 
the normalized GPD H(x,0,t)/q(x) will be constant in 
the limit of x -> 1 |5ITU| . 

In addition to what was mentioned before, GPDs of- 
fer probably for the first time an undisputed way to learn 
about the contribution of quark orbital angular momen- 
tum (OAM) to the nucleon spin. Starting from the QCD 
energy momentum tensor and applying the standard Noether 



At this point it is essential to note that AIq(0) is equal 
to the momentum fraction (x) q ' 9 while -B|q 9 (0) is the sec- 
ond moment of the GPD E, B q ' 9 (0) = / dxxE q > 9 {x, 0, 0). 
The vanishing of the so called anomalous gravitomagnetic 
moment -Bfo~ S (0) has also been derived in the light-front 
formalism 12 . Given the two Eqs.©, the nucleon spin 
sum rule Eq. seems to be a trivial identity. This is only 
true as long as we do not decompose into quark and gluon 
as well as spin and orbital degrees of freedom. A standard 
way of rewriting the spin sum rule is jllj 

\ = \M + L« + J 9 = \ (A w (0) + B 20 (0)) 9+9 , (7) 

where we did not further split the gluon piece into spin 
and orbital angular momentum, since it is not clear how 
to define L 9 properly on an operator level. Assuming that 
a commonly accepted and reasonable definition for the 
gluon spin Ag exists, we may however just set L 9 = J 9 — 
^ffdl- This ambiguity does not affect the following dis- 
cussions, since our current nucleon structure lattice com- 
putations are restricted to quark observables anyway. From 
Eq. Q we find that quark OAM in the nucleon can be com- 
puted once the quark momentum fraction {x} q , the quark 
spin Aq and the second moment of the GPD E q (x,0,0) 
are known. 

Let us now comment briefly on the relation of the 
magnetic moment [i respectively F 2 (0) to quark OAM. It 
has been pointed out frequently J2j that a non-vanishing 
-^2(0) requires necessarily the change of orbital quantum 
numbers Al = ±1 in the corresponding wave functions. 
It is tempting to assume that for this reason a large pro- 
ton (flavor combination 2/3u — l/3<i) magnetic moment 
H p = 1 + F|(0) = 1 + J dxEP(x, 0,0) implies a substan- 
tial quark OAM contribution to the nucleon spin. This is 
not the case. The relevant flavor combination for the spin 
sum rule is of course u + d, and assuming iso-spin sym- 
metry, the experimental value of the iso-sing let F£ +d (0) 
is rather small and negative w —0.36 due to the fact that 
the separate flavors are of almost equal magnitude but op- 
posite sign. In addition, the value for the second moment 
B 2 l Q d (0) is probably even smaller by a factor w 2 ... 4 or 
more ^2- Hence the total quark OAM in the sum rule 
Eq.lJTJ) receives little contribution from the GPD E |15| . 
Even for the separate flavors u and d, it is by no means 
clear that relatively large values of lead to large L u ' d , 
since there may be cancellations of (x) u,d , Aq u ' d and B^ 1 
inL q = l/2{(x)+B 20 ~Aq). 

In summary, the large proton magnetic moment is not 
related to the size of the quark OAM contribution L q to 
the nucleon spin and therefore of no help with regard to 
the " nucleon spin crisis" . This does not contradict the fact 
that we need a change of orbital quantum numbers in the 
wave function to get a non-zero F 2 (0). 
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2 Outline of Lattice Computation of GPDs 

Since it is not feasible to evaluate matrix elements of bilo- 
cal light-cone operators directly on the lattice, the first 
step is to transform the LHS in Eas. l|ll2|l to Mellin space 
by integrating over x, i.e. J dxx n ~ x . The resulting matrix 
elements of towers of local operators are in turn parame- 
terized in terms of generalized form factors (GFFs), here 
generically denoted by A(t), B(t) . . ., 

(P', A'| O r \P, A) = (P', A'| ^(0)riD^- ■ ■ iD^^{0) \P, A) 
= U(P', A') A{t) + bp~» n B(t) + ---)U(P, A), (8) 

where we do not explicitly show the subtraction of traces. 
The explicit parameterization for the tower of vector oper- 
ators, i~=7 /1 , is presented in JS] in terms of the (2n + 1) 
independent GFFs A ni (t), B ni {t) and C n o(t). In the ax- 
ial vector case, the corresponding expression in terms of 
the (2 [I J + 2) GFFs A m (t) and B m (t) is shown in \TQ. 
Recently, parameterizations for the tensor GPDs became 
available as well |17U18j . allowing for a first-time compu- 
tation of these objects in lattice QCD. The GFFs which 
parameterize the matrix elements depend only on the in- 
variant momentum transfer squared t. Therefore, having 
knowledge only of the GFFs, the directional dependence 
via £ = — n ■ A/2 seems to be lost and there is no imme- 
diate way to tell where the information on the transverse 
or longitudinal structure is hidden. It turns out, however, 
that the moments of the GPDs and the GFFs are related 
by comparatively simple polynomial relations in £ which 
restore the full directional dependence J5| • For the vector 
operator we have for example that 



H n (£,t) = / dxx n - x H{x,i,t) 



n-l 
i— O.cvcn 



(-20 1 A ni (A 2 



+ {-2£) n C n0 (A 2 )\ 

n-l 

E n (U)= £ {-20 l B m {A 2 ) 



(-20" C n0 (A 2 )\ 



(9) 



showing that for a purely transverse momentum trans- 
fer, £ = 0,t = —A\, the complete nucleon structure 
concerning H is represented by the set of GFFs A n fi(t), 
n = l...oo. The well known form factors F\ = A\q, 
qa — A\q and gr = Atiq give in this respect just trans- 
verse distributions of quarks in the nucleon, integrated 
over quark longitudinal momentum x. 

On the lattice side, we have to study and eventually 
compute nucleon two- and three-point correlation func- 
tions. The three-point function is defined by 

C pt (r, P', P) = J2 ^ (^nk, i y )Or(r)JV i (r« tc> P)) , 

3,k 

(10) 

where r is a (spin-)projection matrix and the operators 
N and N create respectively destroy states with the quan- 
tum numbers of the nucleon, and Or stands in general for 
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Fig. 2. The lowest three moments of H(x, £ — 0, t), normalized 
to A d j° lc (0) = 1 



the tower of local operators we would like to investigate, 
inserted at a certain timeslice r. To see the relation of 
Co to the parameterization Eq. (JSJ respectively the GFFs 
most clearly, let us rewrite En. I|l()fl using complete sets of 
states and the time evolution operator, 

C 3 pt (r,P',P) = AT(p' ) P) e -^)(----)-E(P')(^-T) 
x ( p ', A 'l Or \P A) U(P, X)fU(P', A') + . . . (11) 

A, A' 

and N(P', P) denotes some normalization factors. The 
dots in (|llfl stay for excited states with energies E 1 > E 
which are exponentially suppressed as long as r — r src > > 
1/E', t sic — t >> 1/E'. Inserting Eq.JSJ, we can carry out 
the sums over polarizations and get 



'(r,P',P) a Tr 



r(p' + m) 



>A(t) 



bp 



>B(t) + •••)(/' + m) 



(12) 



where we dropped for simplicity all pre-factors but showed 
the dependence of the three-point function on the indices 
(/til . . . fi n ). The trace in Ea. (|12J) can be explicitly evalu- 
ated, while the normalization factor and the exponentials 
in Ea. 1)11(1 will be canceled out by constructing an appro- 
priate r- independent ratio R (see e.g. [211] ) of the two- 
and three-point functions. On the lattice side, the ratio 
R is numerically evaluated and then equated with the 
corresponding sum of GFFs times P and P'-dependent 
calculable pre-factors, coming from the trace in Ea. (|12|) . 
For a given moment n, this is done simultaneously for all 
contributing index- combinations (jii ...(j, n ) and all lat- 
tice momenta P, P' corresponding to the same value of t. 
Following this, we end up with an in general overdeter- 
mined set of equations from which we finally extract the 
GFFs. Many details related to the actual calculations can 
be found in |2()j . 



3 Numerical results 

We use for our computations O(200) SESAM configura- 
tions with rather heavy (unimproved) Wilson fermions at 
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m„=897 MeV u-d 
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Fig. 3. The lowest three moments of H(x, £ = 0, t), normalized 
to A d j° lc (0) = 1 

a coupling of /3 = 5.6. Although our results are based on 
a dynamical QCD calculation (including fermion loops), 
we do not take into account the so-called disconnected 
contributions which are numerically much harder to eval- 
uate. We therefore prefer to show the iso-vector channel 
u — d, for which the disconnected pieces, due to iso-spin 
symmetry, would cancel out anyway. The lattice size is 
16 3 x 32, and the physical lattice extension is L ks 1.6fm. 
The results shown below corresponds to a pion mass of 
w 897MeV (k = 0.1560), equivalent to a pion Comp- 
ton wavelength of w 1/7 of the lattice extension. The 
charge radius of the proton turns out to be » 0.5fm for 
these lattice parameters. A selection of our results on the 
^dependence of the GFFs is shown in Figs.ipEjl [2"Tll2"2"]. 
In all cases we used as a lowest order approximation a 
dipole to fit the lattice numbers. The values for the corre- 
sponding dipole masses are 

m D (Ai... 3i o) = 1.47±.03, 2.10±.08, 3.86±.49 GeV 

m D (X..3, ) = 1.85±.03, 2.22±.06, 2.23±.09 GeV 

m D (Si... 3 ,o) = 1.16±.02, 1.66±.02, 2.28±.09 GeV 

(13) 

It is astonishing how strong the slope in t decreases, go- 
ing from the first to the third moment of H, Fig.®. 
But this is what we generically expected from the first- 
principles discussion below Eq.J3J: The higher the mo- 
ment n, the more weight lies on values of x close to one. 
In the extreme case, linXn-^oo A n o(t) / A n o(0) — 1, inde- 
pendent of t. Since the average value of x of the distri- 
bution q u - d {x) = H u - d {x,0,0) is w 0.22 and changes for 
x 2 q u ~ d (x) (third moment) just to « 0.4 [53], the effect ap- 
pears to set in quite promptly |24|. As Fig.(j2Jl indicates, 
the picture looks different for the helicity-dependent GPD 
H . We still see a change in the dipole mass going from A w 
to A30 , but it is much less pronounced than in the helicity 
averaged case. If this picture stands up to further studies, 
one could draw significant conclusions about the helicity 
and transverse momentum structure of the nucleon wave 
function from this observation. The physical interpreta- 
tion of E and the corresponding GFFs B„i concerning the 
transverse nucleon structure is more involved Still, it 
is interesting to note that from the dipole fit, we get a 
noticeable change in the slope in t for the lowest three 
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Fig. 4. The lowest three moments of E(x, £ = 0, t), normalized 

to ^r c (o) = i 

moments as well, Fig.Q. The lattice results for the helic- 
ity dependent analogue, E, suffer at least for the higher 
moments from strong fluctuations and are not shown in 
this work. 

Concerning quark OAM, we have plotted in Fig.JS} 
£?20 an d B d separately. Note that in this case the discon- 
nected pieces have been neglected, leading to a systematic 
error of unknown size. The dipole fits give 

B «,di P oio (0) = o 29± .03, B 2 d dipol ° - -0.38±.02 (14) 

This is exactly what we expected from the discussion at 
the end of section Q). The up and down contributions 
are similar in size but opposite in sign, giving a small 
and negative iso-singlet B%£ d (0) = -0.09 ± 0.04, almost 
compatible with zero. Using the results for the momentum 
fraction and the quark spin contribution |23j . we get for 
quark OAM, L q = \/2{{x) + B 20 - Aq), 

L u = -0.088±0.019, L d = 0.036±0.013 . (15) 

Both flavors separately give a rather small contribution 
of the order of 17% (7%) for u-quarks (d-quarks) to the 
nucleon spin, indeed due to cancellations in quark momen- 
tum fraction, spin and B2o- Adding u and d contributions 
gives a very small and negative total OAM. We would like 
to remind the reader that our calculation correspond to 
a situation where the pion weighs 897 MeV. The corre- 
sponding result for the iso-vector magnetic moment, for 
example, is ~ 20% and the quark momentum fraction 
rs 40% above the experimental value at this pion mass, 
so that the numbers have to be taken with appropriate 
caution. A detailed comparison with experimental results 
calls for lattice computations at much lower pion masses 
and an appropriate chiral extrapolation (work in progress, 
see e.g. |2"5]). 

Lattice calculations of the two lowest moments of the 
helicity independent GPDs in quenched QCD have also 
been presented in Calculations in the quenched ap- 
proximation of quark angular momentum including dis- 
connected contributions have been reported in [2Tll2"8"| . 

4 Summary 

The impact of the observation that GPDs are directly re- 
lated to coordinate space distributions of partons in the 
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Fig. 5. B20 and B20 with dipole fits 



nucleon is profound and even led to the notion of " femto- 
photography" and "tomographic images" of the nucleon 
29 ,30]. Certainly a lot of effort still has to be put into 
the study of the impact parameter dependence of GPDs 
in order to get a detailed picture of the 3-dimensional nu- 
cleon structure. Acknowledging that it is up to now hard 
to come by experimental results on GPDs via deeply vir- 
tual Compton scattering and similar setups, interesting 
analyses making use of experimental data on form fac- 
tors and structure functions to describe GPDs have been 
published recently fR] ( see also As we have shown, 

however, lattice QCD as a first principle technique plays 
already now an important role in unraveling main quali- 
tative features of GPDs and thereby of the nucleon wave 
function. As soon as the infrastructure allows for exten- 
sive calculations using chiral fermions, we expect that the 
methods presented here will give an essential contribution 
to the quantitative understanding of the nucleon structure 
in terms of GPDs. Promising attempts in this direction are 
under way in the framework of a hybrid calculation using 
Astaq sea quarks and domain- wall valence quarks • 
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